SYMBOLIC MODELS FOR NONLINEAR TIME VARYING TIME DELAY SYSTEMS 
VIA ALTERNATING APPROXIMATE BISIMULATION 



GIORDANO POLA 1 , PIERDOMENICO PEPE 1 AND MARIA D. DI BENEDETTO 1 



Abstract. Time-delay systems are an important class of dynamical systems that provide a solid mathemat- 
ical framework to deal with many application domains of interest. In this paper we focus on nonlinear control 
systems with unknown and time-varying delay signals and we propose one approach to the control design of 
such systems, which is based on the construction of symbolic models. Symbolic models are abstract descrip- 
tions of dynamical systems where one symbolic state and one symbolic input correspond to an aggregate of 
states and an aggregate of inputs. We first introduce the notion of incremental input— delay— to— state stability 
and characterize it by means of Lyapunov-Krasovskii functionals. We then derive sufficient conditions for the 
existence of symbolic models that are shown to be alternating approximately bisimilar to the original system. 
Further results are also derived which prove the computability of the proposed symbolic models in a finite 
number of steps. 



keywords: Time-delay systems, symbolic models, alternating approximate bisimulation, incremental input- 
delay— to-state stability, time-varying delays. 

1. Introduction 

Time-delay systems are an important class of dynamical systems that provide a sound mathematical frame- 
work to deal with many application domains of interest, from biology, to chemical, electrical, and mechanical 
engineering, and economics (see e.g. |Nic01j 1. Over the years, several researchers focused on stability, stabi- 
lization, regulation and linearization problems for time-delay systems (see e.g. [HL931 INicOll IKM991 [Ric03 
and the references therein). However, the increasing complexity of current technology asks for control algo- 
rithms that can also deal with a different kind of specifications, such as safety and liveness properties, obstacle 
avoidance, fairness constraints, language and logic specifications, among many others (see e.g. |TP06| ITab09] ) . 
In this paper, we consider nonlinear control systems with unknown and time-varying delay signals and we 
propose a method to deal with such specifications which is based on symbolic models. Symbolic models are 
abstract models where each symbolic state and each symbolic label represent an aggregation of continuous 
states and an aggregation of input and delay signals in the original model. Since these symbolic models are 
of the same nature as the models used in computer science to describe software and hardware, they provide a 
unified language to study problems of control where software and hardware interact with the physical world. 
Moreover, the use of symbolic models allows one to leverage the rich literature developed in the computer 
science community, as for example supervisory control |RW87| and algorithmic game theory [AVW03j . for 
control design of purely continuous processes. Many researchers have recently faced the problem of identifying 
classes of control systems admitting symbolic models. For example, controllable linear control systems and 
incrementally stable nonlinear control systems were shown in |TP06j and respectively in }PGT081 IPT09| , to 
admit symbolic models. Symbolic models for multi-affine systems have been proposed in BH06 and for non- 
linear switched systems in |GPT10| . Symbolic models for nonlinear systems with known and constant delays 
have been explored in [PP DT10] , In this paper we extend the results of the work in PPDT10 to nonlinear 
control systems with unknown and time-varying delays. This class of time-delay systems arises in modeling 
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many application domains of interest and are also challenging from a theoretical point of view, because of 
their inherent complexity, see e.g. |FS06| and the references therein. We stress that the extension of a control 
theory for systems described by time-delay functional differential equations, from the case of known time- 
delays to the case of unknown time-delays, is in general not straightforward, as well as the extension from the 
case of constant time-delays to the case of time-varying time-delays. As far as the stabilization of general 
systems with time-varying time-delays is concerned, by the use of control Lyapunov functionals, the reader 
can refer to the recent book |KJ11| . The methodology of control Lyapunov functionals, in general, requires 
the knowledge of the time-delays. Most nonlinear control laws available in the literature, mainly dealing with 
the stabilization and the input-output linearization problems, are in general designed under the assumption 
of knowledge of the time-delay (often assumed constant), see e.g. [JanOll IGMPOOj IOWN02| IMMM04) . A 
big effort is being made by researchers in recent years to deal with unknown or time-varying time-delays, 
for instance to manage networking control, where the communication time-delay is often time-varying and 
unknown, while an estimation of its upper and lower bounds can generally be given with a certain confidence 
(see for instance |GHL05] for an adaptive controller in the case of unknown constant delays). It is shown 
in |KJ11[ IKJIOj that suitable triangular time-delay systems can be stabilized by a delay-free controller, thus 
completely solving the stabilization problem for this class of systems, also in the case of unknown time-varying 
time-delays. However, to our knowledge, the control of general nonlinear systems with unknown time-varying 
time-delays is still an open and challenging research problem in many cases, for instance when the above 
mentioned specifications are required. 

The framework considered in this paper requires the results of [PPDT10] to be extended in several direc- 
tions. First, an appropriate generalization of the notion of incremental stability proposed in |PPDT10] is 
needed, in order to quantify the mismatch of state trajectories in the symbolic models and in the time-varying 
time-delay systems. To this aim the notion of incremental input-delay-to-state stability (<5-IDSS) is here 
introduced and characterized in terms of Lyapunov-Krasovskii functionals. The incremental input-delay- 
to-state stability (<5-IDSS) is a novel notion that characterizes the relationship between solutions obtained 
by different time-delay signals, besides by different initial conditions and inputs. Second, the presence of 
the additional time-varying time-delay signal in the differential equation requires appropriate approximating 
schemes in order to incorporate its effects in the symbolic model. We solve this problem by approximating 
this time-delay signal by a (quantized) first-order spline [Sch73], that then can be regarded as an additional 
unknown disturbance label in the symbolic model. To our knowledge, it is the first time that the (unknown) 
time-delay signal is approximated by piece-wise linear functions, for the purpose of building up a suitable 
controller. Generally, for this aim, the only state space is projected into suitable finite dimensional subspaces 
(see |GMP00] and references therein). Third, an appropriate notion of approximate equivalence is needed in 
order to capture the adversarial nature of the control input labels and the aforementioned disturbance labels. 
To this aim, we resort to the notion of alternating approximate bisimulation, recently introduced in |PT09| . 
that guarantees that control strategies synthesized on the symbolic models, based on alternating approximate 
bisimulations, can be readily transferred to the original model, independently of the particular realization of 
the time-varying delay signals. While the work in [PPDT10] shows the existence of symbolic models that are 
approximately bisimilar |GP07j to incrementally input-to -state stable nonlinear control systems with constant 
and known time-delays, the results in the present work show existence of symbolic models that are alternating 
approximately bisimilar [PT09 to incrementally-input-delay-to-state stable nonlinear control systems with 
time-varying and unknown time-delays. Further results are also derived which prove the computability of 
the proposed symbolic models in a finite number of steps. An illustrative example is finally presented which 
shows the control design of a 2-dimensional time-delay system with a synchronization-type specification. A 
preliminary investigation on the existence of symbolic models for time-varying time-delay systems appeared 
in the conference publication [PPplCT. In this paper we present a detailed and mature description of the 
results announced in |PPD10| . including proofs and an example. 

The paper is organized as follows. In Section 2 we introduce the class of time-delay systems under con- 
sideration. In Section 3 we introduce and characterize the notion of incremental-input-delay-to-state sta- 
bility. In Section 4 we recall the notions of approximate (bi-) simulations and of alternating approximate 
(bi-)simulations. Section 5 is devoted to the study of existence of symbolic models and Section 6 provides 
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constructive results of the proposed symbolic models. In Section 7 we present an illustrative example. Section 
8 offers some concluding remarks. For the sake of completeness, a detailed list of the employed notation is 
included in the Appendix (Section |9j. 



2. Time- Varying Time-Delay Systems 

In this paper we consider the following nonlinear time-varying time-delay system: 
, s f x(t) = f{x(t),x(t-A(t)),u(t-r)), i£M+«.e. 

[ ' 1 z(t) = &(*). *e[-A max ,o], 

where: 

• x(t) G R n and x t G C7°([— A max , 0]; R n ) is the state at time t G Up ; we recall that x t {6) = x(t + 6) 
with 9 G [-A max , 0] and t G M^; 

• Co = € C°([— A max , 0]; M n ) is the initial condition; 

• u(t) G K m is the control input at time t G [— r, +oo[ and r G Rq is the constant control input delay; 

• A : Rg" — > [A m j n , A max ] is the unknown time-varying delay with A m ; n , A max G R^; 

• / : R n x R" x R m — > R" is the vector field. 



Assumption (A.O). In this paper we suppose that: 

• Function / satisfies /(0,0, 0) = and it is Lipschitz on bounded sets, i.e. for every bounded set 
K C R" x R" x R m , there exists a constant k£I+ such that 

\\f(xi,yi,ui) - f(x 2 ,y2,u 2 )\\ < k(\\x 1 - x 2 \\ + \\yi - 2/2 1| + - "all), 

for all (#1,2/1,111), (x 2 , V2,u 2 ) G if; 

• ||tt(*)|| < 5(7, Vt G [— r, +oo[ for some known B v G M + ; 

• Delay bounds r, A m i n and A max are known; 

• The initial condition £ G C 1 ([-A laax , Q];R n ) with 

(2-2) su PrehAmax , 0] ||e (r)|| < B%, su Pre[ _ Amax!0] ||f (r)|| < M u 

where B% G M + and Mi elj are known; 

• function A G V is unknown and 2? is the set of all continuously differentiable functions A : Rq — > 
[A min , A max ] with ||A(i)|| < d min , where d min G [0, 1[ is known. 



We denote by U the set Bb v (0) and by U the class of measurable control input signals from [— r, +00 [ to U. 
Moreover the symbols x(t, £0, u, A) and Xt(£o, u, A) den ote the solution at time t G in R™ and respectively 
in C°([— A max , 0]; M. n ) of the time-delay system in (2.1 1 with initial condition £0, input u EU and time-delay 
signal AeP. In the sequel we refer to a time-delay system as in (2.1) and satisfying Assumption (A.O), by 
means of the tuple: 



(2.3) E = (M n ,C°([-A E 

where each entity has been defined before. 



,0};R n ),Zo,U,U,V,f), 



3. Incremental Input-Delay-to-State Stability 



The results of this paper rely upon some stability notions that we introduce and characterize in this section. 
We start by recalling the notion of input -to-state stability. 
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Definition 3.1. Son89l [PJ06| A time-delay system £ is Input -to-State Stable (ISS) (uniformly with respect 
to the time-delay function) if there exist a ICC function Piss and a K function 7/55 such that for any time 
t £ Kg , any initial condition £ € C°([— A max , 0]; K ra ), any input 14 eW, and any time-delay function AeP. 



the solution :rt(£g,u, A) of (2.1 1 exists for any t > and, furthermore, the following inequality holds: 
(3-1) ||a: t (£o,u, A)^ < max{/3 J . ? <j(||£o|| 00 ,*),7/ss(||w|[-r,t-r[|| tx) )}- 

A characterization of the ISS property in terms of Lyapunov-Krasovskii functionals can be found in |PJ06j . 
In the sequel, we suppose that: 



Assumption (A.l). The time-delay system in (2.3) is ISS. 

Under the above assumption any solution x(t, £q> it, A) (the initial condition is here considered as part of the 
solution) belongs to B Bx (0) for any time t > -A max , with B x = ma,x{f3 ISS (B^- 7 0),-f I ss{B u )}, where f3 IS s 
and are the ICC and respectively, the K. function appearing in ( |3.1| ). We denote by X the set Bb x (0) 
and by X the set C°(\— A max , 0]; X). Under Assumption (A.l) there is no loss of generality in replacing W 1 



and C ([— A max ,0];M") in (2.3), by X and X respectively, resulting in: 



(3.2) X=(X,X,Z ,U,U,'D,f). 

We can now introduce the notion of incremental input-delay-to-state stability, that adapts the notion of 
incremental input-to-state-stability, introduced in |Ang02| for the class of nonlinear control systems, to the 
class of time-varying time-delay systems considered in this paper. 

Definition 3.2. A time-varying delay system £ satisfying Assumptions (A.0) and (A.l), is incrementally 
Input-Delay -to -State Stable (J-IDSS) if there exist a ICC function (3 and K, functions ju> Id such that for 
any time t £ K^j", any initial conditions £i,£2j any inputs u\,u 2 and any time-delay functions Ai,A2, the 



corresponding solutions Xt(£i, U\, Ai) and Xt(£,2, U2, A2) of (2.1) exist for any t > and, furthermore, the 
following inequality holds: 



(3-3) 1U{\\{U1 -U 2 )\[-r,t-r[\\ 00 ) 

+7d(||(A 1 -A 2 )| [0 , 4[ || oo )}. 

The above notion extends the one proposed in [PPDT10 for time-delay systems with known and constant 
delays, to time-delay systems with unknown and time-varying delays. The main novelty relies on the last 
term in the right-hand side of (3.3). The introduced stability notion regards solutions obtained by different 
time-delay signal realizations, and not only by different initial conditions and inputs. In general the inequality 
in (3.3) is difficult to check directly. We therefore provide hereafter a characterization of 5-IDSS, in terms of 
Lyapunov-Krasovskii functionals. 

Definition 3.3. Given a time-delay system £ = (X, X, £g, U, Li, T>, /), a locally Lipschitz functional 

V:R+x C°([-A maX) 0];R") x C°([-A max , 0];M") -> K+ 

is said to be a 5-IDSS Lyapunov-Krasovskii functional for £ if there exist /Coo functions ai, a 2 and JC functions 
03, p such that: 

(i) for all t € M+ Xl ,x 2 G C°([-A max , 0]; R n ) 

ai(||x 1 (0)-x 2 (0)||) < V(t,xi,x 2 ) < a 2 (M a (x 1 - x 2 )), 
where M a : C°([— A max , 0]; W 1 ) — > Kg is a continuous functional such that 

7 Q (|W0)||)<Af a (z)<7 a (IMU), 
for all x £ C°([— A max , 0]; R™), for some /Coo functions 7 and 7 a ; 
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(ii) by setting L = ma,x{M 1 ,sup xyeX ueU \\f(x,y,u)\\}, ioi all xx, x 2 € X, U\, u 2 EU,de S L(Amax _ Am . n) (0) C 
W 1 for which 

M a (xi - X 2 ) > p{\\{u\ - u 2l d)\\) , 
the following inequality holds, almost everywhere in t € Mj: 

D + V(t,xi,x 2 ,ui,u 2 ,d) < -a 3 (M a (xi - x 2 )), 
where D + V(t,xi,x 2 ,ui,u 2 ,d) is the derivative of the functional V defined as 



V 

D + V(t, X\, x 2 , ui, u 2 , d) — limsup — 

9-»o+ 



(t + e^ 6 ^,^) -v(t,x u x 2 ) 



with 



e,A(t), x _ / Xi(s + 9), if s G [-A max , -9[, 
Xl (S) \ x 1 (0) + (s + 9)f(x 1 (0),x 1 (-A(t)) + d,u 1 ), if sG [-9,0], 

e,A(t)/ s_ j x 2 (s + 9), if S G [-Amax, -6[, 

2 lSj \ s 2 (0) + (« + e)/(ar 2 (0),a: 2 (-A(t)),u 2 ), if* e [-0,0]- 

Before providing a characterization of <5-IDSS in terms of Lyapunov-Krasovskii functionals we need some 
preliminary technical results that we report hereafter. 

Lemma 3.4. Let c e R, p £ N and g : [c, +oo[— > W be a locally absolutely continuous function. Then, for 
an U £i)£2 G [c, +oo[ with ^2 > £i, the following equality holds: 



.9(6) - 



1 dg(0 



o d£ 



dt (6-6). 

«=?i+i(«2-Ci) / 



Lemma 3.5. For any solution x(t,t;o,u, A) attimet> — A max of system (3.2) and for any a, b € [— A max ,+oo[ 
i/ie following inequality holds: 

(3.4) ||a(6 1 f ,u,A)-a;(a,Co,«,A)|| <i|6-o|, 

where we recall L — max{Af 1; sup.,, jj£X uSf7 ||/(x, y, u)||}. 

Proof. By Assumptions (A.O) and (A.l) and the definition of X, any solution lives in X for any time 
t > — A max , and therefore the norm of the derivative of the solution in [— A max ,+oo[ is upper bounded 
by sup X ij/gj y ueU \\f(x, y, u)\\ in R + and by M% in [— A max , 0]. Then, the result holds as a direct application of 
Lemma 13.41 □ 



We now have all the ingredients to present the following result. 

Theorem 3.6. A time-delay system £ is S-IDSS if it admits a 5-IDSS Lyapunov-Krasovskii functional. 

Proof. Let Ai, A2 be two realizations of time-delays, <p\, <p 2 be two initial conditions, u\,u 2 be two inputs. 
Let Xi(t), i=l,2, be the solutions of ( |2.1| ) corresponding to (f> i: u i: A,, respectively. Consider: 

*i(t) = f{xx{t),xx{t - A 2 (i)) + d{t),u x {t - r)), 
d(t) = xi(t - Ai(t)) -xt{t- A 2 (i)), 
x 2 (t) = f(x 2 (t),x 2 (t - A 2 (t)),u 2 (t - r)). 



By Lemma 3.5 it follows that 



(3.5) \\d(t)\\ < L||Ai(t) - A 2 (t)|| < i(A max - A min; 
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By (3.5 1, the proof is over if the existence of the 5-IDSS Lyapvmov-Krasovskii functional guarantees that the 
system described by 

x(t) = f(x(t), x(t - A(t)) + d(t), u(t - r)) 

is such that there exist a ICC function j3, JC functions ju, Id such that, for any realization of the time-delay 
A(t), initial conditions 0j, inputs Ui(t), i = 1,2 and d(t) (with their upper bounds), the following inequality 
holds for the solutions, 



,o r N \\xt(<f>l,Ul,d) - X t ((f>2,U 2 ,0)\\ < ^Bx{P(\\^ 1 -4>2\\oo,t),^u{\\(ul-U2)[- r ,t-r\\ 00 ) 

{ } +^(\\d m \\J}- 

We now show that the same lines of the proof used by Sontag for ISS (see, for instance, |SW95] ). used also 
for time-delay systems in [PJ06L can be used here. Let 0i, <f> 2 be a pair of initial conditions. Let U\, u 2 , d be 
such that || (tii — ^2 , cZ) 1 1 oo = v- The set 

S = {(t^ 1 ,^ 2 )eR + x C°([-A max ,0];R n ) x C°([— A max , 0];R") : V(t,i/? u ih) < <*2 ° p(v)} 

is such that: if (to,x to (<j)i,ui,d),Xt (4> 2 ,u 2 ,0)) € S for some to £ Rjj", then 

(t,x t (<t>i,ui,d),x t (<t>2, u 2 ,0)) € S, 

for all t > to. To prove this, by contradiction suppose that there exists a time i > to such that V{i, Xt{4>\,u\,d), Xt(<fi 2 ,u 2 , 0)) > 
a 2 o p(v). Then there exists a positive real e and a time t\ > to such that 

V(ti,x tl (^i,ui,d),x tl (4> 2 ,u 2 ,0)) > a 2 o p(v) + e. 

Let ti be minimal for this given e. By continuity arguments, it follows that there exists a € M. + such that 
for any t e]ti — a,t± + a[, V(t, xt{4>i, u\, d), Xt(<f> 2 , u 2 , 0)) > a 2 o p(v). By (ii), Theorem 5 in |Pep07b] and 
Theorem 2 in |Pep07a| , it follows that V(t, Xt((f>i, u\, d), Xt(4> 2 , u 2 , 0)) is nonincreasing in ]t\ — a,t\]. Thus t\ 
is not minimal and a contradiction holds. Now two cases can occur: 1) (0, 01,02) belongs to S; 2) (0,01,0a) 
does not belong to S. In case 1) by the forward invariant property of the set S, the following inequality holds, 
for any t > 0, 

\\x(t, 4>i,ui,d) ~ x(t, 02, u 2 , 0)|| < a^ 1 o a 2 o p(v). 
In case 2) let to & M(| U {+00} be the maximal time (notice that it can be +00) such that, in the interval 
[0, to[) (t, Xt(4>i, Mi, d), Xt{4> 2l u 2 , 0)) does not belong to S. Then, in [0,io[> the inequality in (ii) holds, that, 
by Theorem 2 in |Pep07a| , yields for w(t) — V(t,x t ((j>i,ui,d), x t ((f> 2 ,u 2l 0)), the inequality D + w(t) < — a 3 o 
a 2 1 { w {t)) a - e - From this inequality, by (i), Theorem 5 in |Pep07b| , the following inequality holds, for a suitable 
ICC function (3 (see Lemma 4.4 in |LSW96j ). for t G [0,t o [, 

(3.7) \\x(t, 0i,ui,d) - x(t, 02,^2, 0)|| < a^ 1 o j3(a 2 °7 a (||0i - 02||oo),£)- 

By the result concerning the forward invariant property of the set S, the following inequality holds, Vi > 0: 

, . ||a:(t,0i,ui,d) - x(t, 2 , u 2 , 0)|| < max{a^ o (3(a 2 o 7 a (||0i - 2 ||oo), t), 

-1 ( \1 



one 



Therefore, in both cases 1) and 2), the inequality (3.8 1 holds for any t £ R . From the inequality (3i 
gets: 

\\x t ((f>i,ui,d) - x t (02,u 2 ,O)||oo < max{e~ ( *~ Amax) ||0i - 2 ||oo 
(3-9) +ar 1 °^(«2O7 Q (||0i_-0 2 ||oo), 

max{0, t - A max }), 0/0020 p(v)} 

By causality arguments (sec SW95 ), i.e. the solutions x t ((f>i,Ui,d) and x t (4> 2 , u 2 ,0) are independent of the 
inputs ui, u 2 at time greater than t — r (recall the input delay), and of the input d at time greater than t, it is 



possible to replace v by ||((ui — u 2 )\[_ r t _ r [, cf| ro,t[) lloo - Thus, the inequality in (3.6) is proved, which concludes 
the proof. □ 
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Remark 3.7. Sufficient and necessary conditions for a time-delay system to be ISS, in terms of existence of 
ISS Lyapunov-Krasovskii functionals, can be found in |KPJ08b] and in |QGD10| . For a time-delay system 
to be (5-IDSS it is unclear at this stage whether the existence of a (5-IDSS Lyapunov-Krasovskii functional is 
not only a sufficient condition but also a necessary one. When inputs are not involved (i.e. = U2(t) = 0, 

Vi > 0), and the same realization of time-delay signal is considered (i.e. Ai(i) = ^(t), Vf > 0), the 
notion of (5-IDSS reduces to the well known <5-GAS notion (see |Ang02| ), extended to nonlinear time-delay 
systems. The (5-GAS notion can be also defined as a GAS notion (see |An g02 ) with respect to a suitable 
diagonal set by means of an extended system constructed as follows: the original system is replaced by 
two subsystems, each one described by the same dynamic functions of the original one. To this aim, the 
following result, that can be regarded as an extension of Lemma 2.3 in Ang02|, is useful. Let the product 

l ) be endowed with the norm given, for ip — 



space S = C°([-A 



V>i e C°([-A max ,0]; 

" " 



as 



4> 



,0];R n ) x C°( 
"), %= 1,2, as 
,iP€C°([-A max ,0};R n ) 



-A max ,0];R n ) be endowed with the norm given, for = I e S, 
0||§ = ||^i||oo + Halloo- Let H be the subset (diagonal set) of S defined 
. Then, for any 0i,0 2 € C°([-A mox , 0];R n ), the equality holds 



(3.10) 



||-02 - "01 



inf 

X ec°([-A ma 



01 

i>2 



x 
x 



The proof of (3.10) is as follows. Take any ip 



(3.11) 

On the other hand. 
(3.12) 



inf 

X ec°([-A ma 



,0];K") 



01 







X 
X 



01 
02 

< 



G S. The inequality holds 







01 
01 



< inf ( 

X eC°([-A max ,0];K") 



= ||02 ~ 0l||oo- 
02 - Xlloo + 1101 - X\\oo) 



inf 

X6C°([-A mox ,0];K") 



0- 



X 
X 



From (3.11 1 and (3.12) the equality (3.10) follows. This result is in no contradiction with t he result in Lemma 
2.3 in Ang02| (where a factor y/2/2 is involved). Actually the same result shown in (3.10) holds for the finite 
dimensional case, when the norm in R 2rl = M n x R" is chosen as the sum of the euclidean norms in R™. 

01 



Notice that the right-hand side of (3.10) is a point distance of 



02 



from the diagonal set EL Therefore, 



the (5-GAS of the original system is equivalent to GAS property of the above mentioned extended system, 
with respect to the diagonal set H, also for time-delay systems. By using this result, the converse Lyapunov- 
Krasovskii theorems available in the literature for the GAS property (see |KP J08a) . |K J 11] ). extended to the 
case of stability with respect to the diagonal set, can be used directly to provide converse Lyapunov-Krasovskii 
theorems for the <5-GAS notion of time-delay systems. Whether the results in [KPJ08b could be used also 
for the (5-IDSS, in order to provide converse Lyapunov-Krasovskii theorems, is an interesting and not easy 
topic of future investigations. The main issue concerns dealing with different time-delay signal realizations, 
which to our knowledge has never been considered in the literature so far. 



Quadratic Lyapunov-Krasovskii functionals GKC03] are popularly used for the stability check of time-delay 
systems. They are an important set of functionals that can be useful also for checking the (5-IDSS property, 
as it will be shown in Section 7 through an illustrative example. 



4. Symbolic Models and Approximate Equivalence Notions 



We use the class of alternating transition systems [AHKV98 as abstract mathematical models of time-delay 
systems. 
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Definition 4.1. An (alternating) transition system is a tuple: 

T=(Q,q ,L, » ,0,H), 

consisting of: 

• A set of states Q; 

• An initial state go € Q; 

• A set of labels L = A x B, where: 

— A is the set of control labels; 

— B is the set of disturbance labels; 

• A transition relation ► C Q x L x Q; 

• An output set O; 

• An output function H : Q — > O. 

A transition system T is said to be: 

• metric, if the output set O is equipped with a metric d:OxO-> Rq~ ; 

• countable, if Q and L are countable sets; 

• finite, if Q and L are finite sets. 

We follow standard practice and denote by q °' » ■ p, a transition from q to p labeled by a and b. Transition 

systems capture dynamics through the transition relation. For any states q,p £ Q, the symbol q — p simply 
means that it is possible to evolve or jump from state q to state p under the action labeled by a and b. 
We consider simulation and bisimulation relations [Mil89, Par81 as formal equivalence notions to relate time- 
delay systems and their symbolic models. Intuitively, a bisimulation relation between a pair of transition 
systems T\ and T 2 is a relation between the corresponding state sets explaining how a state trajectory T\ of 
T\ can be transformed into a state trajectory r 2 of T 2 , and vice versa. While typical bisimulation relations 
require that r\ and r 2 have the same output run, i.e. H 1 (ri) — H 2 {r 2 ), the notions of approximate bisimulation 
relations, that we now introduce, relax this condition and require that Hi(r\) is simply close to ^2(^2) where 
closeness is measured with respect to a metric on the output set. A simulation relation is a one-sided version 
of a bisimulation relation. 

Definition 4.2. [GP07] Given two metric transition systems T\ = {Q\,q1,A\ x B\, — — -+■ , 0,B\) and 

T 2 = (Q2, q 2 ,A 2 x B 2 , — iO,H 2 ) with the same output set O and the same metric d and given a precision 

e £ M.q, a relation R C Q 1 x Q 2 is said to be an £-approximate simulation relation from T\ to T 2 if for any 
(qi,q 2 ) E R: 

(i) d(H 1 (q 1 ),H 2 (q 2 ))<e; 

(ii) V(ai,6i) £ Ai x B Y , ^(a 2 ,b 2 ) £ A 2 x B 2 such that q% pi and q 2 p 2 with (pi,p 2 ) £ R. 

Moreover, T± is said to be e-approximately simulated from T 2 , denoted T\ < e T 2 , if there exists an e- 
approximate simulation relation R from T\ to T 2 such that {q®, q 2 ) £ R. Given a precision e £ Kj, a relation 
R C Qi x Q 2 is said to be an e-approximate bisimulation relation between T\ and T 2 if: 

(i) R is an e-approximate simulation relation from T± to T 2 ; 

(ii) is an e-approximate simulation relation from T 2 to T\. 

Moreover, T\ is said to be e-bisimilar to T 2 , denoted T\ = e T 2 , if there exists an e-approximate bisimulation 
relation R between T\ and T 2 such that {q\, q 2 ) £ R. 



The above notions have been further extended in |PT09| to the ones of alternating approximate simulation 
and bisimulation, that appropriately capture the different role played by control labels and disturbance labels. 
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Definition 4.3. Given two metric transition systems Ti = (Qi,q®,Ai x £> l5 — — »- , 0,Hi) and T 2 — 

(Q2,'?2!^2 x B2, — with the same output set O and the same metric d and given a precision 

e G Kg", a relation R C Q 1 x Q2 is said to be an alternating e-approximate (AeA) simulation relation from 
T\ to T 2 if for any {q\,q%) G -R: 

(i) d(/Ti(gi),JT 2 (g 2 )) <e; 

(ii) Vai G Ai 3<22 G A 2 V62 G B2 36i G -Bi such that 51 ^ »- pi and q2 *■ P2 with (pi, 1*2) G i?. 

Moreover, Ti is said to be AeA simulated from T2, denoted Ti T2, if there exists an AeA simulation 
relation R from Ti to T2 such that (q®, q%) G R. Given a precision e G Rq , a relation i? C Qi x Q2 is said to 
be an alternating e-approximate (AeA) bisimulation relation between Ti and T2 if: 

(i) R is an alternating e-approximate simulation relation from Ti to T 2 ; 

(ii) is an alternating e-approximate simulation relation from T2 to Ti. 

Moreover, Ti is said to be AeA bisimilar to T2, denoted Ti =| T2, if there exists an AeA bisimulation relation 
R between Ti and T 2 such that (q°,q°) G R. 

When e = 0, the notion of bisimulation in the above definition coincides with the 2-players version of the 
definition proposed in [AHKV98 . 

5. Existence of Symbolic Models 

5.1. Alternating Approximately Bisimilar Symbolic Models. In this section we propose symbolic mod- 
els that approximate (5-IDSS time-delay systems in the sense of AeA bisimulation. These symbolic models 
are obtained by a first time discretization of the time-delay system and a subsequent quantization of the state 
and input variables, and of the delay signals. 

In many concrete applications controllers are implemented through digital devices and this motivates our 
interest for control inputs that are piecewise constant. In the following we refer to time-delay systems with 
digital controllers as digital time-delay systems. We suppose that the set U of input values contains the origin 
and that the class of control inputs is 

(5.1) U T := {u : [—r, — r + r] — > U \ 3v G U s.t. u = u|r_ ri _ r _|_ T i}, 

where r G K + can be thought of as a sampling time parameter in the digital controller. For further purposes, 
define: 

V T := {A : [0, r] -> [A min , A max ] | 3d G V s.t. A = d| [0 , T] }. 
Given k G N we denote by Uk, T and by T>k, T the class of control inputs and time-delay signals obtained by 
the concatenation of k control inputs in U T and respectively, of k time-delay signals in T> T . Given a digital 
time delay system S define the transition system 

T T (S) -^(Q^q^L,, ,01,/Ti), 

where: 

• Qi = X; 

• q°i - £0; 

• Li = A\ x i?i, where: 

— A\ = {a\ G U T I x T {x, a\, b\) is defined for any x G X and b\ G B{\\ 

— Bi = {&i G T> T I x T (x, ai, 61) is defined for any x G X and ai G j4i}; 

• q P, if x T {q, ai, bx) = p; 

• 0\ = X: 
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• H x =l x . 

Transition system T T (S) can be thought of as a time discretization of S and it is metric when we regard 
0\ = X as being equipped with the metric d(p, q) = \\p — q\\oo- Note that T T (S) is not symbolic, since the set 
of states Qi is a functional space. The construction of symbolic models for digital time-delay systems relies 
upon approximations of the set of reachable states, of the space of time-delay signals and of the space of input 
signals. Let i? r (E) C X be the set of reachable states of £ at times t = 0, r, fcr, i.e. the collection of 
all states x £ X for which there exist k £ N, a control input u £ U^,t an d a time-delay signal A £ T>kri 
so that x — Xk T (£,o,u, A). The sets i? T (£), A\ and B\ are functional spaces and are therefore needed to be 
approximated, in the sense of the following definition. 

Definition 5.1. jPPDTlOj Consider a functional space y C C°{I,Y) with Y C E", I = [h,i 2 ], h,i% £ K, 
i\ < i 2 . A map A : M + — > 2 C is a countable approximation of y if for any desired precision A £ R + : 

(i) -4(A) is a countable set; 

(ii) Vy e J, 3z e .4(A) so that \\y - z\\oo < A; 
(hi) Vz e -4(A), By £ y so that jjy - zjjoo < A. 



A countable approximation Au of U T can be easily obtained by defining for any \jj £ IB 
(5.2) Au(\u) = {u£U T : u(t) = «(-r) G [ET^,* G [-r.-r + r]}, 



where [I^Alt is defined as in (9.1 1. By comparing W r in (5.1) and Ajj{Xu) in ( |5.2| it is readily seen that 
Au{^u) C W T for any Ay £ M + . Since the set U contains the origin, the set Au{^u) is nonempty. Suppose 
that the time-delay system £ is 5-IDSS and let j3, ju and be a K.C and K. functions satisfying the inequality 
in ( |3.3[ ). We can now define a countable transition system that approximates T T (E). Given a time quantization 
r elF, a state space quantization Ax £ R + , an input space quantization Xjj £ R + and a delay quantization 
\d £ consider the transition system: 

(5-3) T t , Ax ,a^a d (£) := (Q2,q° 2 ,L 2 , ,0 2 ,H 2 ), 



where: 



• q% £Q 2 such that ||£ - q 2 \\oo < Ax; 

• L 2 = A 2 x U 2 , where 4 2 = Au(\jj) and B 2 = Ad(A d ); 

• 9 — ^ P, ii 

||p- x T (g,a 2 ,6 2 )|| 00 < max{/3(Ax,r),7 D (Ai))} + Ax; 

• 2 = 4"; 

• H 2 = i : Q 2 ^ 2 - 

We can now give the following result: 

Theorem 5.2. Consider a digital time-delay system £ and an?/ desired precision e £ R + . Suppose that £ is 
S-IDSS and choose r € M + so t/iai /3(e, r) < e. Moreover suppose that there exist countable approximations 
Ax and Ad of R t (E) and T> T , respectively. Then, for any Ax £ R + , A/y £ R + and Ad £ R + satisfying the 
following inequality: 

(5.4) max{/?(e, r), 7t/(Aa) + 7d(Xd)} + max{/3(Ax, r), 7d(A jD )} + Ax < e, 

transition systems T T \ X \ U \ D (E) andT T (E) are As A bisimilar, i.e. T Tt \ Xj \ Ut \ D (E) T T (E). 



Before giving the proof of the above result we point out that: 

Lemma 5.3. For any given precision e £ K + , there exists a choice of quantization parameters r, Ax, \u> A^ £ 
R + so that the inequality in (5.4) holds. 
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Proof. Pick Ax < e /3- Since 7^7 and 7^ are K, functions, there exists a choice of Au and A^ so that "fu(^u) < 
e/6 and 7,d(A_d) < e/6. Since /3 is a ICC function there exists r so that /8(e,r) < e/3. By construction, since 
Ax < e, then max{/3(Ax, r), 7_d(Ad)} < max{/3(e,T),7i/(Ai/) + 7d(Ad)}, from which: 

max{/3(e,T),7( 7 (A C /) + Jd{Ad)} + max{/3(Ax, r), 7d(Ad)} + Ax < 
2max{^(e,r),7[/(A[ / ) + 7d(A d )} + Ax < 2max{e/3,e/6 + e/6} + e/3 = e, 



which concludes the proof. 



□ 



We can now give the proof of Theorem |5.2 



Proof. Consider the relation R C Q 1 x Q 2 defined by (x, 9) £ R if and only if ||i? 1 (x) — #2(9)11 00 < e - We now 
show that R is a AeA simulation relation from T T (E) to T Tt \ x \ Ut \ D (?l). Consider any (x,q) G R. Condition 
(i) in Definition 4.3 is satisfied by the definition of R. Let us now show that condition (ii) in Definition 4.3 
holds. Consider any a\ G A\. By definition of A2 there exists 02 G A2 so that: 

(5.5) ||a x -CI2IL < At/. 

Consider any 62 G B 2 . By construction of T T ^\ x x u x D {^) there exists 61 G Bi so that: 

(5.6) ||foi - &2II00 < Azj. 

By definition of Q2 there exists x € i? r (E) so that: 

(5.7) ||2-9||oo < Ax- 

Set y = x r (x, 02, 61) G i? r (£). By definition of Q2 there exists p € Q2 so that: 
(5-8) ||tf-p||oo < A x . 

Set z = x T (g, a2, 62)- Note that since a2 G A2 C £/ r , function z is well defined. By the <5-IDSS assumption 
and the inequalities in (5.7) and (5.8), the following chain of inequalities holds: 



(5.9) 



\z-p\\<x> = \\z - y + y -p\\oo 

< \\z- y\\oo + \\y-p\\oo 

< max{/3(||g - x\\ 00 ,T),j u (\\a2 - a 2 ||oo) + 7.d(||&i - &2II00)} + Ax 

< max{/3(Ax,r),7 C/ (0)+7 D (A I) )} + Ax. 



By the above inequality it is clear that q — p in T Ti \ Xi \ Vj \ d {Ti). Consider x ai ' l y. Since E is 5-IDSS 



and by (5.4), (5.5), (5.6) and (5.9), the following chain of inequalities holds: 



\y-p\\oo < \\y - z\\oo + \\z - p\\oo 

< max{/3(||x - 9||oo,T),7t/(ll a i ~ a 2 \\) + 7z>(||i>i - &2IU)} 
+ max{/3(Ax,r),7 £) (AD)} + Ax 

< max{^(e, r),7t/(At/) + Jd{>>d)} + max{/3(Ax, r), 7d(Ad)} + Ax < e. 



Hence (y,p) G i? and condition (ii) in Definition 4.3 holds. By the inequality in (5.4) and the definition of q®, 
1 1 £0 — <?2 1 — Ax < £ and hence, transition system T T (S) is AeA simulated by T T} x x ,\ U} x D (^)- By using similar 
arguments it is possible to show that R^ 1 is an AeA simulation relation from T Tj x x ,Xu,*d (^) t° ^V(S). Hence, 
the result follows. □ 
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5.2. Spline-based countable approximation of functional spaces. The result presented in the previous 
section assumes existence of countable approximations of functional spaces of time-delay systems. By following 
PPDT10J in this section we present an approach to approximate these functional spaces which is based on 
spline analysis |Sch73j . Spline based approximation schemes have been extensively used in the literature of 
time—delay systems (see e.g. |GMP00j and the references therein). Let us consider the space y C C°(I,Y) 
with Y C K™, I = [f i, £2], ii,i2 G K and i\ < i^. Given N <E N consider the following functions (see |Sch73j ): 

1 — (t — h)/h, t€[ii,ii + h], 
0, otherwise, 



s {t) 



(5.10) 



l — i + (t — ii)/h, te [h + (i — h + ih], 

s l (t)={ l+i-\t-ii)/h, te [ii+i/i,ii + (i + l)/i], 

0, otherwise, 
1 = 1,2,. ..,7V; 



SjV+l(*) 



0, otherwise, 



where h = {ii — ii)/(N + 1). Functions Sj called splines, are used to approximate 3^. Given any JV e N. 
(9, M € R+ let b«Q 

(5.11) A(N, 6, M) := h 2 M/8 + (N + 2)6, 

with h = (%2 — ii)/(N + 1). Function A will be shown to be an upper bound to the error associated with the 
approximation scheme that we propose. It is readily seen that for any A e K + and any M e K + there always 
exist N e N and 9 e K+ so that A(N, 9, M) < A. Let Nx,m and 9 x ,m be such that A(N\ tM ,8\,M,M) < A. 
For any A e K + and M e R + , define the operator ^a,M : y Y) that associates to any function y e y 

the function: 

Nx.M + l 

(5.12) i>\,M{y)(t) := £ te-T, 

j=0 

where ^ 6 [^]2e A M and — y(«i + < 0\.m, for any i = 0, 1, N\,m + 1. Note that the operator ip\ t M 
is not uniquely defined. For any given M € K + and any given precision A € K + define: 

(5.13) ^,m(A) := Va,m(30- 

Lemma 5.4. [PPDTlOj Suppose that y C PC 2 {I;Y) and there exists M e R+ so that ||£> 2 i/||oo < M for 



any y e J 7 - TTien ^lj;,Af as defined in ,5.13), is a countable approximation of y . 

The above result is useful for approximating the functional space T> T of time-delay signals, as shown hereafter: 
Proposition 5.5. Consider the functional space T> T and Mu e K + . Suppose that: 

(A.2) V T = ?C 2 ([0,r]; [A min , A max ]); 
(A.3) for any A e V T , \\D 2 A\\ <X> < M D . 

Then the set Ad defined for any \r> G M + by: 

(5.14) A D (\ D )=i>x D ,M D (V T ), 



with ip\ Dt M D as i n (5-12), is a countable approximation ofV T . 



We now proceed with a further step towards the construction of countable approximations of i? T (S). Consider 
a digital time-delay system X = {X,X,^,U,U T ,V T ,f) and suppose that: 

(A.4) £ is 5-IDSS; 



^The real M is a parameter associated with y and its role will become clear in the subsequent developments. 
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(A.5) Function / is Frechet differentiable in W 1 x E" x R m ; 

(A. 6) The Frechet differential J(x, y, u) of / is continuous and bounded on bounded subsets of R" x R" x M™. 
Under the above assumptions the following bounds are well defined: 

(5.15) Bj = s\ip^ >VtU - )eXxXxU \\J{x,y, u)\\, M x = {2B X + B v )(l + d min )KBj, 

where k is the Lipschitz constant of function / in the bounded set X x X x U and || J(x, y, u)\\ denotes the 
norm of the operator J(x, y, u) : K" x W 1 x M m — > R n . We can now give the following technical lemma that 
is instrumental to prove the main result of this section. 

Lemma 5.6. Consider a digital time-delay system £ satisfying assumptions (A.O), (A.l), (A.^-6), and 

(A. 7) the following conditions: 

(A.7.1) Co G C°([-A max ,0];X)nPC 2 ([-A max ,0];]R™); 

(A.7.2) H^olL <M X ; 

(A.7.3) P I ss(B x ,t)+ 1isS (B u ) < B x ; 

(A.7.4) t > 2A max . 

Then, for any x kr € i? r (S) with k G N the following hold: 

• x kT G C°(hA max ,0];A)nPC 2 ([-A max ,0];]R"); 

• H^ferlloo < B° x ; 

• \\D 2 x kT l <M X . 



Proof. It is sufficient to show that x T satisfies the same properties of £o> i- e - conditions (A.7.1), (A.7.2) hold 
with £g replaced by x T and ||^ r ||oo < B x . First note that the function t — > x(t), t G [0,r], is uniformly 
continuous in the (compact) set [0,t]. By Assumption (A.7.4), it follows that x T+ g G C 1 ([— A max , 0]; E"), 
8 G] — A max ,0[, i.e. the derivative i T +e G C°([— A max , 0]; R n ). By taking into account the Lipschitz constant 
k (computed on the bounded set X x X x U) of function /, the bounds B x and Bjj, the following chain of 
inequalities holds: 

ll-ir+elloo = su Pae[-A max ,o] \\f(x(T + e + a),x(T + e + a- A(r + 6» + a)), 
u(t + 9 + a - r))|| 
< KSup ae[ _ Amax 0] (||a;(r + 6 + a)\\ + \\x(t + 9 + a - A(r + 6 + a))||+ 
||it(T + + a-r)||) 



< k(2Bjc + £c), fle]-A 



,0 



As far as the second derivative is concerned, the following equality holds, for g] — A max ,0[, 



<je 2 



J(x(r + 8),x(t + 8-A(t + 6)),u(t + 8- r))< 
i(r + 0) 

f(r + 6» - A(r + 6>))(1 - A(r + 8)) 


Since the Frechet differential J is continuous, since x T+ e G C°([— A max , 0]; M n ) and from the differentiability 
assumption on the state time-delay, it follows that x T G PC 2 ([— A max , 0]; M n ). Moreover, by taking into 
account the bound £>,/ on the Frechet differential J, the bound on the derivative £ T +g in (5.16) and the 
definition of M x in (5.15), we obtain ||Z3 2 x T || 00 < M x . Thus x T satisfies conditions (A.7.1) and (A.7.2). 
Finally by condition (A.7.3) it is readily seen that H^Hoo < B x . □ 



The above result shows that under assumptions (A. 0-1), (A. 4-6) and (A. 7. 3-4), the regularity properties of 
the initial state £o in (A. 7. 1-2) propagate to the whole set of reachable states, or in other words, that time- 
delay systems are invariant with respect to the properties in (A. 7. 1-2). This is a key result that allows us 
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to underline sufficient conditions for the existence of a countable approximation of i? r (£), as formally stated 
hereafter. 

Theorem 5.7. Consider a digital time-delay system S, satisfying assumptions (A. 0-7). Then the set Ax 
defined for any Xx G R + by: 

(5.16) Ax(Xx) = rpx x ,M x (Rr(X)), 



with ipx x ,M x as i n (5.12), is a countable approximation of i? T (S). 



The proof of the above result is a direct consequence of Lemmas |5.4| and |5.6| and it is therefore omitted. 

5.3. Main Result. We now have all the ingredients to define a symbolic model for digital time-delay systems. 
Given r G M + , 8x,0d, Xjj,Md G K + and Nx, Np G N, consider the transition system 

(5-17) T T,(N x ,6 x ),\ u ,(N D ,e D ){?') ■= r T ,A x ,A[/ 

where T T ,x x ,Xu,X D (£) is defined in (5.3) with Xx = A(N X ,6 X ,M X ) and Ad = A(Nd, 0d, Md). It is readily 
seen that Assumptions (A. 0-7) guarantees that transition system T T f^ X! g x \x tr ,m D ,e D )0^) m (5.17) is symbolic. 



We can now present the main result of this paper. 

Theorem 5.8. Consider a digital time-delay system £ and any desired precision e G K + . Given Md G R + , 
suppose that assumptions (A. 0-7) are satisfied. Moreover let t,9 x ,9d, Xu G R + and N x ,Nd G N satisfy the 
following inequality 

max{/3(e, r), 7c/ (A l /) + ^jd(A(Nd, #d, 

max{/3( A( A x , X , M x ), r), lD (A(N D , D , M D ))} + A(N X , 6 X , M x ) < s, 



(5.18) 



with A as in (5.11) and M x as in (5.15). Then transition systems T T (£) and T T /N Xt g x \x v JND,6D)C^) are 
AsA bisimilar, i.e. T T . (Nxfix)Au ^ ND j D )(Y,) =f T T (E). 

Proof. The map Au is a countable approximation of U, by Proposition |5.5| the map Ad is a countable 



approximation of D T and by Theorem 5.7 the map Ax is a countable approximation of i? T (E). Choose 
Ax G R + , Xrj G R+ and Ay G R + satisfying the inequality in ( 5\4| . There exist 0x G M + and Ax G N so 
th at Ax = A(N X ,6 X ,M X ), 9 D G R + and A^ G N so that X D = A(N D ,9 D ,M D ), that satisfy the inequality 
in (5.18). Finally the result holds as a direct application of Theorem 5.2 □ 



6. Construction of Symbolic Models 



The construction of the symbolic model T , r (x x ,0 x ).A (7 ,(A' D ,e o )(£) in (5.17) requires the preliminary computa- 
tion of countable approximations Au and Ad oiU T and T> T , respectively. While the computation of Ajj is 
straightforward, the computation of Ad is not so, because Ad is defined as the image through the operator 
ipX D ,M D °f the set T> T that is composed by an infinite and uncountable number of functions. In this section 
we present some results that are weaker than the one in Theorem |5.8| for which in turn, symbolic models can 
be effectively constructed. The main idea is to define a suitable symbolic model T*, Nx g ^ ^ g ^{Ti) that 
can be effectively computed and that approximates the symbolic model T T (x x ,e x ),A [/ ,(Af D ,e £ ,)(S), in the sense 
of alternating 0-approximate simulation relation. Given 

Tt,(N x M x ),X u ,(N d ,8 d ){^) = (Q2,q2,L 2 , —r* ,0 2 ,H 2 ), 



with L 2 
(6.1) 



A 2 x B 2 , define the symbolic model: 



T- 



r,(N x ,9 x ),X u ,(N D ,0 



d) (S) :=(Q 2 ,q° 2 ,L* 2 , ,0 2 ,H 2 ), 



where L 2 = A 2 x B 2 and B 2 is the collection of all functions b(t) 



jji G [D]g D . Note that by definition B 2 C B 2 from which, the following result holds. 



,0], with 
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Proposition 6.1. 

T t,(N x ,9 x ),\u ,(N D ,6»d)( S ) — 0* T T ,(N x ,6x),>>u,(ND,6n)(^) -0 T t,(N x ,6 x ) ,\ v ,(N d ,6 d ) 



The proof of the above result is a straightforward consequence of the definition of the symbolic models involved 
and it is therefore omitted. We now have all the ingredients to present the following result. 

Theorem 6.2. Consider a digital time-delay system E and any desired precision e € K + . Given Mrj £ K + , 
suppose that assumptions (A. 0-7) are satisfied. Moreover let r, Ox,&d, € M + and Nx,Nd € N satisfy the 
inequality in \5.18 ), with A as in (5.11) and Mx as in (5.15). Then: 



Proof. (Proof of T* 
(6.2) 



6.1 



i.Vx,«,A„( Wl ,,fo)( E ) T ^ E ))' ^ ^position 

^(JVx.ex)^^,^,^)^) -o t ^V,(Af x ,ex),A !7 ,(Af D ,e £ ,)(E). 



By Theorem 5.8 and since AeA bisimulation implies AsA simulation, one gets: 
( 6 - 3 ) t t,(Nx ,e x )M ,(N D ,0d)(£) r<£* ?V(E). 



Hence, by combining (6.2 1 and (6.3) and by a straightforward generalization of Proposition 2 in |GP07j . the 
result follows. The proof of the approximate inclusion T T (E) ;< £ T* , Nx fl x Xu , Nd q d \{£) can be given along 
the lines of the proof of Theorem 5.8 and is therefore omitted. □ 



From I AHKV98 , PT09J , the above result guarantees that control strategies synthesized on the symbolic model 
-^r (N x e x ) \u (N D D )(^) can ^ e reao -^y transferred to the original system E, independently of the particular 
realization of the time-varying delay signal A. The above result is weaker than Theorem |5.8| in the sense 
that it does not guarantee existence of alternating approximate bisimulation between the time-delay system 
E and the corresponding symbolic model T* {Nx ^ x)Au ^ Nd ^ d )( 



5.8 



does). The motivation 



(E) (as Theorem 

in the introduction of the symbolic model T* , Nx e ^ x ^ e D )(^) ^ s that ^ can be effectively computed, as 
discussed hereafter. It is easy to see that the set B 2 coincides with the co-domain of the operator ipo D: M D ( scc 
(5.12) and (5.14)) and it is composed by a finite number of functions. Hence, the set B 2 can be computed in a 
finite number of steps from which, the symbolic model in (6.1 ) can be effectively constructed. The construction 



of the proposed symbolic models can be easily derived by adapting Algorithm 1 in |PPDT10] for symbolic 
models of time-delay systems with constant delays to symbolic models of time-delay systems with time- 
varying delays. However, the adaptation of Algorithm 1 in [PPDT10] to this framework is not efficient from 
the computational complexity point of view because it generally leads to large symbolic models and extensive 
time of computation that are not needed when solving many (symbolic) control design problems. A more 
efficient approach would construct the symbolic controller without constructing the whole symbolic model of 
the time-delay system. Useful insights in this direction are reported in |PBD12j . concerning the integrated 
symbolic control design of nonlinear systems. We do not report in the paper technical results generalizing the 
ones of |PBD12] to time-delay systems for lack of space. We instead illustrate in the next section, through a 
simple example, the computational complexity gain obtained by following this approach. 



7. AN ILLUSTRATIVE EXAMPLE 



In this section we illustrate the results presented in this paper by means of a simple example. Consider the 
following nonlinear time-delay system: 

x x {t) = -8xi (t) + tanh(x 2 (t - A(t))) 
x 2 (t) = -9x 2 (t) + sin(xi(t - A(t))) + cos(x 2 (t))u(t - r), 



(7.1) 
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GIORDANO POLA, PIERDOMENICO PEPE AND MARIA D. DI BENEDETTO 



where t G 



A(t) G 



with A, 



10" 



10- 



10, d n 



0.2, and tanh(x) 



e»+e~ x ^ 0r an y ^ S K. We address a symbolic control design problem where the state of E is requested 
to reach some regions of the state space within some prescribed times. More specifically, we consider the 
following synchronization specification: starting from the origin, remain in the positive orthant for all times; 
reach the set X\ = [0.01, co[x [0.1, oo[ in no more than 4s, stay in the set X\ for at least 4s, reach the set 
X 2 = [0.08,015] x [0.08,0.15] and finally remain in X 2 for at least 12s. These requirements arise for example 
in multi-agent systems with shared resources in which the use of a given resource is needed to be synchronized 
among the agents. More complex specifications can be also considered as logics-based specification, fairness 
constraints, and etc. (see e.g. |TP06j ). In order to solve the considered control design problem we first need 
to check 5-IDSS of S. For the system to satisfy Assumption (A.l), consider the quadratic functional (see 
Remark 3.9 in |PJ06j . concerning the linear increasing kernel in the second integral term) defined, for any 



t G 



^0 ! 



G C°([-A max ,0];M 2 ), 



as: 



Vr SS (t,0)=^(O) + ^(Q) 



(7.2) 



+ 



-A„ 



-A(t) 

-TAT r Q (T + A max ) 



>>\(t) + 4>\{T))dT 



A, 



A, 



4> 2 (r)dr, 



where r n , r\ are positive reals, with r n > r&. By applying Theorem 3.1 in pJ06| and appropriately choosing 
the parameters ro, ta, we obtain that Viss is an ISS Lyapunov-Krasovskii functional for system (7.1|. In 



particular, for ta = 0.2, r$ = 0.3, the inequality in (3.1| is fulfilled with functions: 

Pi SS {u,t) = 2.8920e- 1 - 0870t w, jiss(w) = 0.9592 uj u,t€ Rq . 

We now proceed with a further step and consider Assumption (A. 4). Consider the quadratic functional defined, 
for any t G K+ Mfa G C([— A max , 0]; M 2 ), as: 



V(t, 0!, 0a) = (<M0) - M0)f(M0) - M0)) 

+2 f (Mr) ~ Mr)f(Mr) ~ Mr))dr 
J-A(t) 

''" / -tat r (r + A 

max J 



A maJ 



A, ; 



(t) - Mr)V (Mr) ~ Mr))dr, 



(7.3) 



where Tq, ta are positive reals, with tq > ta- By appropriately choosing the parameters ro, r a, we obtain that 



V is a <5-IDSS Lyapunov-Krasovskii functional for system (7.1) and hence, by Theorem 3.6 the time-delay 



system E is 5-IDSS. In particular, for ro = 0.3, ta = 0.2, the inequality in (3.3 1 is fulfilled with functions 



P{u,t) = 4.3580e- 1 - 087C V 7[7 (w) = 13.5647 ui, ~/ d {uj) = 194.1666w, w,(eR+. 

Let be B° x = 0.5, B v = 0.3, Mi = 0.1, and M D = 0.001. Consequently, we obtain Bj = 27.9, k = 
9.3, M x = 993.8845. For a precision e = 0.12, we can choose r = 2, X x = 0.02, X v = 5 • 10~ 4 , and 



Xr, = 1.4 -10 4 so that the inequality in (5.4) is satisfied. By the definition of function A in (5.11) and 



since X x = A(N X ,9 X ,M X ) and X D = A(N D ,9 D ,M D ), one can choose N x = 0, 6 X = 0.04, N D = 1 



and 8d = 6 • 10~ 6 . This choice of quantization parameters satisfies the conditions in (A. 7. 3) and (5.18) 
By generalizing the algorithms proposed in [PPDTIO] . the estimated time to construct the symbolic model 
(N x o x ) Ac/ (N D 8 D )C^) 1S about 203, 215s. Since the expected time of computation is rather high, in the 
sequel we adapt the algorithms proposed in [PBD12] towards the efficient symbolic control design of nonlinear 
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time-delay systems. The symbolic control strategy obtained is reported hereafter, where (m, n-i) 



(n+,n+) 



(7.4) 





)) Ul 


^ {{n+8 x 


n+8 x ) 




Ox)) 


(0,0) 


186 


(4,30) 


-396 


(3,22) 


248 


(4,31) 


-562 


(3,20) 


-268 


(3,24) 


-546 


(2,20) 


-484 


(3,21) 


388 


(4,33) 


234 


(4,31) 


-220 


(3,25) 


542 


(4,35) 


-560 


(3,19) 


-74 


(3,27) 


-142 


(3,26). 





The running time needed for solving the given symbolic control design problem is 7, 692s by using a laptop 
with CPU Intel Core 2 Duo T5500 @ 1.66GHz. Figure [T] shows the evolution of the state variables of E with 
the unknown time-delay signal: 



(7.5) 



A(t) 



2 2 

It is readily seen that the synchronization specification is indeed satisfied 



sin(O.OH), t G 




0.1 



0.05 




10 



15 



20 



25 



30 



Figure 1. Evolution of the state variables of E with initial condition £o = 0, time-varying 



delay signal in (7.5|, and the control strategy reported in (7.4 1 



Conclusion 



In this paper we studied existence of symbolic models for nonlinear control systems with time-varying and 
unknown time-delay signals. We first introduced the notion of <5-IDSS and provided a characterization of 
this property, in terms of Lyapunov-Krasovskii functionals. We then showed that 5-IDSS time-delay systems 
admit symbolic models that are alternating approximately bisimilar to the original systems with a precision 
that can be rendered as small as desired. We finally presented a symbolic model that is shown to be an 
alternating approximate simulation of the original system and that can be computed in a finite number of 
steps under a boundedness assumption on the sets of states, inputs and delays of the system. 
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9. Appendix: Notation 

The symbols N, Z, R, R + and Rq denote the sets of natural, integer, real, positive and nonnegative real 
numbers, respectively. Given a vector x G R n the z— th element of x is denoted by xf, furthermore ||a;|| denotes 
the infinity norm of x; we recall that ||a;|| := max{|a;i|, \x 2 \, \x n \}, where \xi\ is the absolute value of Xi. For 
a given x G R n and a given s G R + , the symbol B s (x) denotes the set {y G R n : \\y — x\\ < s}. For any iCM" 
and 9 G R + define 

(9.1) [A] e := {a e A \ai = kid, fc* G Z,i = 1, ...,n}. 

Given a measurable and locally essentially bounded function / : Rq — > R™, the (essential) supremum norm of 
/ is denoted by ||/||oo; we recall that ||/||oo : = (ess) sup{||/(i)||, t > 0}. For a given time r G R + , define f T 
so that f T (t) = f(t), for any t G [0, r[, and f(t) — elsewhere; / is said to be locally essentially bounded if 
for any r G R + , f T is essentially bounded. A continuous function 7 : Rq~ — > Rq is said to belong to class JC 
if it is strictly increasing and 7(0) = 0; 7 is said to belong to class JC^ if 7 G K, and 7(7*) — > 00 as r — > 00. 
A continuous function j3 : Rq x Rq — > Rq is said to belong to class 1CC if for each fixed s, the map /?(r, s) 
belongs to class K, with respect to r and, for each fixed r, the map f3(r, s) is decreasing with respect to s and 
(3{r,s) -> as s -> 00. Given k,n G N with n > 1 and / = [a, b] C R, a, b G R, a < b let C k (I;R n ) be the 
space of functions / : I — >• R™ that are continuously differentiable k times. Given k > 1, let PC k (I;R n ) be 
the space of C k ~ 1 (I; R n ) functions /:/—>■ R n whose fc-th derivative exists except in a finite number of reals, 
and it is bounded, i.e. there exist 70,71, ■■■,'Js G R + with a = 70 < 71 < ... < 7 S = b so that D k f is defined 
on each open interval (7i,7»+i), i = 0, 1, s - 1 and max i=0 ,i,... lS -i sup te(7t 7 . +i) \\D k /(t)||oo < 00. For any 
continuous function x(s), defined on —A < s < a, a > 0, and any fixed t, < t < a, the standard symbol xt 
will denote the clement of C°([-A, 0];R") defined by x t (6) = x(t + 9), -A < 6 < 0. The identity map on a 
set A is denoted by 1a- Given two sets A and B, if A is a subset of B we denote by %a ■ A B or simply by 
1 the natural inclusion map taking any a G A to i{a) = a G B. Given a function / : A — > B the symbol f(A) 
denotes the image of A through /, i.e. f(A) := {b G B : 3a G A s.t. b = /(a)}. If C C A, the symbol f\c 
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denotes the restriction of function / to C, i.e. f\c( c ) — f( c ) f° r an y c £ C. Given a relation R C Ax B, R 1 
denotes the inverse relation of R, i.e. R^ 1 := {(b, a) G B x A : (a, b) E A x B}. 
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